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Abstract 



The presently most accurate prediction for the W-boson mass in the Standard 
Model is obtained by combining the complete two-loop result with the known higher- 
order QCD and electroweak corrections. The numerical impact of the different 
contributions is analysed in detail. A simple parametrisation of the full result is 
presented, which approximates the full result for M\y to better than 0.5 MeV for 
10 GeV < Mh < 1 TeV if the other parameters are varied within their combined 2a 
region around their experimental central values. The different sources of remaining 
theoretical uncertainties are investigated. Their effect on the prediction of M\y is 
estimated to be about 4 MeV for M H < 300 GeV. 
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The relation between the W-boson mass, Mw, the Z-boson mass, M%, the fine structure 



constant, a, and the Fermi constant, G^, 



Mi 1 - — ^ = 1 + Ar 



is of central importance for precision tests of the electroweak theory. Accordingly, a lot 
of effort has been devoted over more than two decades to accurately predict the quantity 
Ar, which summarises the radiative corrections, within the Standard Model (SM) and 
extensions of it. 

The one-loop result pQ can be written as 

Ar^ = Aa - ^Ap + Ar rcm (M H ), (2) 
s w 

where c w = M w /M|, s w = 1 — c w . It involves large fermionic contributions from the shift 
in the fine structure constant due to light fermions, Aa oc logm/, and from the leading 
contribution to the p parameter, Ap. The latter is quadratically dependent on the top- 
quark mass, m t , as a consequence of the large mass splitting in the isospin doublet jSJ. The 
remainder part, Ar rem , contains in particular the dependence on the Higgs-boson mass, 
M H . Higher-order QCD corrections to Ar are known at 0(aa s ) [3] and 0(aa 2 ) jHIH]- 

Recently the full electroweak two-loop result for Ar has been completed. It consists of 
the fermionic contribution [q T| T || 5] , which involves diagrams with one or two closed fermion 
loops, and the purely bosonic two- loop contribution 0. 

Beyond two-loop order the results for the pure fermion-loop corrections (i.e. contribu- 
tions containing n fermion loops at ra-loop order) are known up to four-loop order [TU1 . 
They contain in particular the leading contributions in Aa and Ap. Most recently re- 
sults for the leading three-loop contributions of O(G^m^) and 0{G 2 ^a s mf) have been 
obtained for arbitrary values of Mh (by means of expansions around Mh = m t and for 
Mh 3> m t ) [TT] . generalising a previous result which was obtained in the limit Mh = 

Eq. ((TJ) is usually employed for predicting the W-boson mass, 




1 + Ar(M w , M z , Mr, m t , . . . ) }, (3) 



which is done by an iterative procedure since Ar itself depends on Mw- Comparison 
of the prediction for Mw within the SM with the experimental value allows to obtain 
indirect constraints on the Higgs-boson mass. These constraints are affected both by the 
experimental error of Mw and by the uncertainty of the theory prediction. The current 
experimental error of the W-boson mass is SM^ P = 34 MeV The accuracy in the 
measurement of the W-boson mass is expected to improve to about 5M w :p ' Tev/ ' LHC = 
15 MeV jT3] from the measurements at Runll of the Tevatron and the LHC, and to about 
5Mw P ' = 7 MeV at a future Linear Collider (LC) running at the WW threshold |15j . 
The uncertainty of the theory prediction is caused by the experimental errors of the input 
parameters, e.g. m t , and by the uncertainty from unknown higher-order corrections. In 
the global SM fit to all data ^H] the highest sensitivity to Mh arises from the predictions 
for Mw and the effective weak mixing angle at the Z-boson resonance, sin 2 9 e g. 
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In the present paper we combine the various pieces being relevant for the prediction of 
Mw into a common result and analyse the numerical impact of the different contributions. 
Since in particular the electroweak two-loop result is very lengthy and involves numerical 
integrations of two-loop scalar integrals, it is not possible to present the full result in a 
compact analytic form. We therefore provide a simple parametrisation of the full result 
which is easy to implement and should be accurate enough for practical applications. We 
discuss the sources of the remaining theoretical uncertainties and obtain an estimate for 
the uncertainty from unknown higher-order corrections. 

We incorporate the following contributions into the result for Ar, 

^ r = Ar(«) + Ar^ aQ ^ + Ar^ aa '^ + Ar| a ^ + Ar^ + Ar ( - G ^ asm ^ + Ar^ m *\ (4) 

where Ar^ is the one-loop result, eq. (J2J), Ar(° Qs ) and Ar( aa ^ are the two-loop Pj and 
three-loop jUEj QCD corrections, and Arj^ 00 E] and Ar£J are the fermionic 
and purely bosonic electroweak two-loop corrections, respectively. The contributions 
Ar^ G l asm ^ and Ar^ G i m ^ have been obtained from the leading three-loop contributions to 
Ap given in Ref. 

We have not included the pure fermion-loop contributions at three-loop and four- 
loop order obtained in Ref. |10j . because their contribution turned out to be small as a 
consequence of accidental numerical cancellations, with a net effect of only about 1 MeV 
in Mw (using the real-pole definition of the gauge-boson masses). Since the result given 
in Ref. ^U] contains the leading contributions involving powers of Aa and Ap beyond 
two-loop order, we do not make use of resummations of Aa and Ap as it was often done 
in the literature in the past (see e.g. Refs. J7|). Accordingly, the quantity Ar appears 
in eq. (j3J) in fully expanded form. This means, for instance, that we do not include the 
£>(a 3 ) term 3(Aa) 2 Arj^, which can be inferred from the electric charge renormalisation. 
It affects the prediction for Mw by about 1.5 MeV. This shift is however expected to 
partially cancel with the corresponding contributions proportional to (Aa)(Ap)Ar^ and 
(Ap) 2 Ar[^ in an analogous way as for the pure fermion-loop contributions. 

In Table ^ the numerical values of the different contributions to Ar are given for 
M\v = 80.426 GeV [T^]. The other input parameters that we use in this paper are 

mt = 174.3 GeV, m b = 4.7 GeV, M z = 91.1875 GeV, T z = 2.4952 GeV, 
a" 1 = 137.03599976, Aa = 0.05907, a s (M z ) = 0.119, 

Gfj, = 1.16637 x 10~ 5 GeV" 2 , (5) 

where Aa = Aa lcpt + Aa^, and Aai ep t = 0.0314977 [TH]. For Aaj^ we use the value 
given in Ref. [H], Ao£ = 0.027572 ± 0.000359. The total width of the Z boson, T z , 
appears as an input parameter since the experimental value of M z in eq. (J3J), corresponding 
to a Breit-Wigner parametrisation with running width, needs to be transformed in our 
calculation into the mass parameter defined according to the real part of the complex 
pole, which corresponds to a Breit-Wigner parametrisation with a constant decay width, 
see Ref. [7]. It is understood that Mw in this paper always refers to the conventional 
definition according to a Breit-Wigner parametrisation with running width. The change 
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Table 1: The numerical values (xlO 4 ) of the different contributions to Ar specified in 
eq. (JTJ) are given for different values of Mh and Mw = 80.426 GeV (the W and Z masses 
have been transformed so as to correspond to the real part of the complex pole). The 
other input parameters are listed in eq. ©. 



of parametrisations is achieved with the one loop QCD corrected value of the W-boson 
width as described in Ref. [7j. 

Tabled shows that the two- loop QCD correction, Ar^ aas \ and the fermionic elec- 
troweak two-loop correction, Arj^ are °f similar size. They both amount to about 10% 
of the one-loop contribution, Ar^ a \ entering with the same sign. The most important 
correction beyond these contributions is the three- loop QCD correction, Ar^ aas \ which 
leads to a shift in Mw of about —11 MeV. For large values of Mh also the contribution 
/\ r (^a s m*) becomes sizable. The purely bosonic two-loop contribution, Ar^ s \ and the 
leading electroweak three-loop correction, Ar'^" 1 '', give rise to shifts in Mw which are 

significantly smaller than the experimental error envisaged for a future Linear Collider, 

5M e*P,LC = 7 MeV 

Since Ar is evaluated in Tabled for a fixed value of Mw, the contributions Ar(° Qs ) 
and Ar < - QQs * ) are Mn-independent. In the iterative procedure for evaluating Mw according 
to eq. (j3J), on the other hand, also these contributions become M H -dependent through the 
M H -dependence of the inserted M w value. 

The result for Mw based on eqs. 0, (JU) can be approximated by the following simple 
parametrisation (see Ref. [20] for an earlier parametrisation of Mw), 

M w = M w -cidH-c 2 dH 2 + C3dH 4 + C4(dh-l)-c 5 da + c 6 dt-c 7 dt 2 

— c 8 dH dt + c 9 dh dt — c w da s + c\\ dZ, (6) 

where 

dH-tof M " V dh=C M » V, dt=( )'-!, 
V 100 GeV J VlOOGeVy ' V 174.3 GeV/ 

dZ= M \ 1, d«= da, = ^M-l, (7) 

91.1875 GeV 0.05907 ' 0.119 ' v ; 
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Table 2: Shifts in Mw caused by varying Mh by 100 GeV and the other input parameters 
by la around their experimental central values [T^j. The first column shows the full result 
for Mw, while the second column is based on the simple parametrisation of eqs. 
The shifts 5M\y are relative to the value Mw — 80.3799 GeV which is the result for 
Mh = 100 GeV and the central values of the other input parameters as specified in 
eq. p. 



and the coefficients M w , ci, . . . , en take the following values 

M w = 80.3799 GeV, a = 0.05429 GeV, c 2 = 0.008939 GeV, 

c 3 = 0.0000890 GeV, c 4 = 0.000161 GeV, c 5 = 1.070 GeV, 

c 6 = 0.5256 GeV, c 7 = 0.0678 GeV, c 8 = 0.00179 GeV, 

c 9 = 0.0000659 GeV, c 10 = 0.0737 GeV, c n = 114.9 GeV. (8) 

The parametrisation given in eqs. (H3)~© approximates the full result for Mw to better 
than 0.5 MeV over the whole range of 10 GeV < Mh < 1 TeV if all other experimental 
input values vary within their combined 2a region around their central values given in 
eq. ©. 

In Table 121 the full result for M w and the parametrisation of eqs. (jHJ)~@ are com- 
pared with each other. The table shows the shifts in Mw (relative to the value Mw = 
80.3799 GeV, which is the result for Mh = 100 GeV and the central values of the other 
input parameters as specified in eq. ©) induced by varying Mh by 100 GeV and the other 
input parameters by la around their experimental central values |13j . In the example of 
Table 121 where only one parameter has been varied in each row and all others have been 
kept at their central values, the maximum deviation between the full result for Mw and 
the parametrisation of eqs. ©-© is below 0.1 MeV. 

The parametrisation of eqs. ©-© yields a good approximation of the full result for 
Mw even for values of Mh much smaller than the experimental 95% C.L. lower bound 
on the Higgs-boson mass, Mh = 114.4 GeV [2Tj . If one restricts to the region Mh > 
100 GeV, a slight readjustment of the coefficients in eq. (jSJ) yields an even more accurate 
parametrisation of the full result. If eqs. (JHJ), (JZj) ar e used together with the following 
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values of the coefficients, 



M w = 80.3800 GeV, c x = 0.05253 GeV, c 2 = 0.010345 GeV, 

c 3 = 0.001021 GeV, c 4 = -0.000070 GeV, c 5 = 1.077 GeV, 

c 6 = 0.5270 GeV, c 7 = 0.0698 GeV, c 8 = 0.004055 GeV, 

c 9 = 0.000110 GeV, c 10 = 0.0716 GeV, c u = 115.0 GeV, (9) 

the full result for Mw is approximated to better than 0.2 MeV over the range of 100 GeV < 
Mh < 1 TeV if all other experimental input values vary within their combined 2a region 
around their central values given in eq. (J7|). 

From Table 121 one can read off the parametric theoretical uncertainties in the pre- 
diction for Mw being caused by the experimental errors of the input parameters. The 
dominant parametric uncertainty at present (besides the dependence on Mh) is induced 
by the experimental error of the top-quark mass. It is almost as large as the current ex- 
perimental error of the W-boson mass, SM^ P = 34 MeV j!3j . The uncertainty caused by 
the experimental error of m t will remain the dominant source of theoretical uncertainty 
in the prediction for Mw even at the LHC, where the error on m t will be reduced to 
5m t = 1-2 GeV [22] • A further improvement of the parametric uncertainty of Mw will 
require the precise measurement of m t at a future Linear Collider [23 , where an accuracy 
of about 5m t = 0.1 GeV will be achievable |15j . 

We now turn to the second source of theoretical uncertainties in the prediction for Mw, 
namely the uncertainties from unknown higher-order corrections. Different approaches 
have been used in the literature for estimating the possible size of uncalculated higher- 
order corrections |2"Hl25|l27)l lTj. The "traditional Blue Band method" is based on the fact 
that the results of calculations employing different renormalisation schemes or different 
prescriptions for including non-leading contributions in resummed or expanded form dif- 
fer from each other by higher-order corrections. The deviations between the results of 
different codes in which the same corrections have been organised in a somewhat dif- 
ferent way are used in this method as a measure for the size of unknown higher-order 
corrections [21] • In applying this method it is not easy to quantify how big the variety of 
different "options" and different codes should be in order to obtain a reasonable estimate 
of the higher-order uncertainties. As the method cannot account for genuine effects of 
irreducible higher-order corrections, it may lead to an underestimate of the theoretical 
uncertainties if at an uncalculated order a new source of potentially large corrections 
appears, e.g. a certain enhancement factor. 

In Ref. j2E] a different prescription has been proposed, in which for each type of 
unknown corrections the relevant enhancement factors are identified and the remaining 
coefficient arising from the actual loop integrals is set to unity. In Ref. |7] higher-order 
QCD corrections have been estimated in two different ways, from the renormalisation scale 
dependence (in particular taking into account the effect of switching from the on-shell to 
the MS definition of the top-quark mass) and from assuming that, for instance, the ratio 
of the 0(a 2 a s ) and 0(a 2 ) corrections is of the same size as the ratio of the 0(aa s ) and 
0(a) corrections. 

Several of the corrections whose possible size had been estimated in Refs. [23112^ 17] 
have meanwhile been calculated [HUH], and it turned out that the estimates agree rea- 
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sonably well with the actual size of the corrections. This adds confidence to applying the 
same kind of methods also for an estimate of the remaining higher-order uncertainties. 

There are three sources of remaining uncertainties in the prediction for from 
unknown higher-order corrections: 

• The corrections at 0(a 2 a s ) beyond the known contribution of 0{G 2 L a s mf): 

The numerical effect of the 0{G 2 l a s ml ) correction was found to be up to 5 MeV in 
M w for a light Higgs-boson mass, M H < 300 GeV [TT]. This contribution represents 
the leading term in an expansion for asymptotically large values of m t . In the 
calculation of the electroweak two-loop corrections it was found that the formally 
next-to-leading order term of 0{G 2 ^m 2 M^ j ) has approximately the same numerical 
effect as the formally leading term of 0(G 2 l mf) [27J. It can therefore be expected 
that also the formally next-to-leading order term of 0(G 2 ^a s m 2 M%) may be of similar 
size as the leading 0{G 2 L a s mf) term. We therefore assign an uncertainty of about 
3 MeV to the remaining theoretical uncertainties at 0(a 2 a s ) (for M H < 300 GeV). 

• The unknown electroweak three-loop corrections: 

The numerical effect of the 0{G 3 l m^) contribution was found to be small jTT], shift- 
ing Mw by less than 0.3 MeV for Mh < 300 GeV. This shift is significantly smaller 
than the estimate in Ref. [26J. The pure fermion-loop corrections at three- loop or- 
der were found in Ref. [10) to shift Mw by about 1 MeV, which however involved 
an accidental numerical cancellation. It thus doesn't seem to be justified to as- 
sume that all other electroweak three-loop corrections are completely negligible. In 
Ref. [7] it was pointed out that reparametrising the W-boson width, which enters 
the prediction for M w at the two-loop level, by G M instead of a shifts the prediction 
for Mw by about 1 MeV, which is formally an effect of 0(a 3 ). In order to take into 
account uncertainties of this kind (see also the discussion below eq. (JH)) we assign 
an uncertainty of 1-2 MeV to the unknown corrections at 0(a 3 ). 

• The four-loop QCD corrections of 0(aa 3 ): 

The possible effect of the leading term of 0(G^alm 2 ) was estimated in Ref. [2H] to 
be about 1.3 MeV. Employing the known results at lower order of a s and assuming 
a geometric progression yields a slightly larger result. We thus assign an uncertainty 
of 1-2 MeV for the 0(aa 3 ) corrections. 

Adding the above estimates for the different kinds of unknown higher-order corrections in 
quadrature, we find as estimate of the remaining theoretical uncertainties from unknown 
higher-order corrections 

5M$ eo « 4 MeV. (10) 

This estimate holds for a relatively light Higgs boson, Mh ^$ 300 GeV. For a heavy Higgs 
boson, i.e. Mh close to the TeV scale, the remaining theoretical uncertainty is significantly 
larger. 

In Fig. ^we have updated the comparison between the theory prediction for Mw within 
the SM and the experimental value, using the currently most accurate theory prediction 
based on eqs. (j2J), (jlj) and the most up-to-date experimental data [13J. For the theoretical 
uncertainty the estimate of eq. (jlOj) and the parametric uncertainties corresponding to la 
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Figure 1: Prediction for Mw in the SM as a function of Mh for m t = 174.3 ± 5.1 GeV. 
The current experimental value, M^ p = 80.426 ± 0.034 GeV and the experimental 
95% C.L. lower bound on the Higgs-boson mass, M H = 114.4 GeV are also indicated. 



variations of the input parameters (see Table EJ) have been used. As discussed above, at 
present the theoretical uncertainty is dominated by the effect of the experimental error 
of the top-quark mass. 

Fig. ^ confirms the well-known preference for a light Higgs-boson mass within the 
SM. If the 95% exclusion bound from the direct search for the SM Higgs is taken into 
account [21], the la bands corresponding to the theory prediction and the experimental 
result for Mw show only a marginal overlap. 

In summary, we have presented the currently most accurate prediction for Mw in 
the Standard Model. We have discussed the relative importance of the complete one- 
loop and two-loop contributions as well as the known corrections beyond two-loop order. 
We have summarised the present status of the theoretical uncertainties of M w from the 
experimental errors of the input parameters, and we have obtained an estimate for the 
remaining theoretical uncertainties from unknown higher-order corrections. In the region 
of Higgs-mass values preferred by the electroweak precision data, M H < 300 GeV, the 
uncertainty from unknown higher-order corrections amounts to about 4 MeV. This is 
much smaller than the present experimental error of Mw and even below the envisaged 
future experimental error at the next generation of colliders. Having reached this level of 
theoretical precision of Mw is important, however, for the precision test of the electroweak 
theory, in particular in view of the fact that Mw can be used as an input for calculating 
the effective weak mixing angle at the Z resonance, sin 2 9 e s. 

We have furthermore presented a simple parametrisation of the full result containing 
all relevant corrections, which should be sufficiently accurate for practical applications. 
It approximates the full result for Mw to better than 0.5 MeV over the whole range of 
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10 GeV < Mh < 1 TeV if all other experimental input values vary within their combined 
2a region around their experimental central values. In view of the experimental exclusion 
bound on the Higgs-boson mass of Mh > 114.4 GeV it will normally be sufficient to 
restrict to the smaller range of 100 GeV < Mh < 1 TeV. For this case we provide a 
simple parametrisation which approximates the full result for M\y even within 0.2 MeV. 
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